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Abstract 

We introduce a "renormalized entanglement entropy" which is intrinsically UV finite and is most 
sensitive to the degrees of freedom at the scale of the size R of the entangled region. We illustrated 
the power of this construction by showing that the qualitative behavior of the entanglement entropy 
for a non-Fermi liquid can be obtained by simple dimensional analysis. We argue that the functional 
dependence of the "renormalized entanglement entropy" on R can be interpreted as describing the 
renormalization group flow of the entanglement entropy with distance scale. The corresponding 
quantity for a spherical region in the vacuum, has some particularly interesting properties. For 
a conformal field theory, it reduces to the previously proposed central charge in all dimensions, 
and for a general quantum field theory, it interpolates between the central charges of the UV 
and IR fixed points as R is varied from zero to infinity. We conjecture that in three (spacetime) 
dimensions, it is always non-negative and monotonic, and provides a measure of the number of 
degrees of freedom of a system at scale R. In four dimensions, however, we find examples in which 
it is neither monotonic nor non-negative. 
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I. INTRODUCTION 



Quantum entanglement has been seen to play an increasingly important role in our un- 
derstanding and characterization of many-body physics (see e.g. [1, 2]). The entanglement 
entropy for a spatial region provides an important set of observables to probe such quantum 
correlations. 

In spacetime dimensions higher than two, however, the entanglement entropy for a spatial 
region is dominated by contributions from non- universal, cutofF-scale physics [3, 4]. This 
implies that for a region characterized by a size R, the entanglement entropy is sensitive 
to the physics from scale R all the way down to the cutoff scale 6, no matter how large R 
is. As a result the entanglement entropy is ill-defined in the continuum limit. The common 
practice is to subtract the UV divergent part by hand, a procedure which is not unique and 
often ambiguous, in particular in systems with more than one scales. Even with the UV 
divergent part removed, the resulting expression could still depend sensitively on physics 
at scales much smaller than the size R of the entangled region. As a result, in the limit 
of taking R to infinity, one often does not recover the expected behavior of the IR fixed 
point (see for example the case of a free massive scalar in Sec. V). 

Such a situation is clearly awkward both operationally and conceptually. We should be 
able to probe and characterize quantum entanglement at a given macroscopic scale without 
worrying about physics at much shorter distance scales. 

In this paper we show that there is a simple fix of the problem.^ Consider a quantum 
field theory on ]R^'°'~^ which is renormalizable non-perturbatively, i.e., equipped with a well- 

^ See also [5] for a discussion based on free theories. 
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defined UV fixed point. Suppose S^^\R) is the entanglement entropy in the vacuum across 
some smooth entangling surface S characterized by a scalable size R.^ We introduce the 
following function 

^(.)(^^ ^ I (^^ - 1) (^il - 3) ■ ■ ■ {R^i -id- 2)) S(^\R) d odd 

1 JI^M (^^ - 2) ■ ■ ■ {R^ -id- 2)) S(^\R) d even 

In Sec. II we show that it has the following properties: 

1. It is UV finite in the continuum limit (i.e. when the short-distance cutoff is taken to 
zero). 

2. For a CFT it is given by a /^-independent constant . 

3. For a renormalizable quantum field theory, it interpolates between the values s^^'^^^ 
and s^^'^^'' of the UV and IR fixed points as R is increased from zero to infinity. 

4. It is most sensitive to degrees of freedom at scale R. 

The differential operator in (1.1) plays the role of stripping from S^^\R) of short-distance 
correlations. The stripping includes also finite subtractions and is /^-dependent; it gets of 
rid of not only the UV divergences, but also contributions from degrees of freedom at scales 
much smaller than R. Sji^\R) can be also be used at a finite temperature or finite density 
where it is again UV finite in the continuum limit. In the small R limit it reduces to the 
vacuum behavior while for large R we expect it to go over to the thermal entropy. 

(R) may be considered as the "universal part" of the original entanglement entropy, a 
part which can be defined intrinsically in the continuum limit. Below we will sometimes refer 
to it as the "renormalized entanglement entropy," although this name is clearly not perfect. 
We believe such a construction gives a powerful tool for understanding entanglement of a 
many-body system. As an illustration, in Sec. Ill we show that the entanglement entropy 
of a non- Fermi liquid also has a logarithmic enhancement just as that for a Fermi liquid by 
a simple dimensional analysis. We also predict the behavior of the entanglement entropy 
from higher co-dimensional Fermi surfaces. 



^ In this paper we will always consider S to be a closed connected surface. Also note that not all closed 
surfaces have a scalable size. In Sec. II and Appendix A we make this more precise. 
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In the rest of the paper we focus on the behavior of (R) in the vacuum, studying 
its possible connections to renormahzation group flow (RG) and the number of degrees 
of freedom along the flow. Items (2)-(4) above, especially (4), indicate that S^\r) can 
be interpreted as characterizing entanglement correlations at scale R. Thus in the contin- 
uum limit as we vary R from zero to infinity, S'^^\R) can be interpreted as describing the 
renormahzation group (RG) flow of the "renormalized entanglement entropy" from short 
to large distances. In contrast to the usual discussion of RG using some auxiliary mass or 
length scale, here we have the flow of a physical observable with real physical distances. Its 
derivative 

(1.2) 

dR ^ ' 

can then be interpreted as the "rate" of the flow. With the usual intuition that RG flow 

leads to a loss of short-distance degrees of freedom, it is natural to wonder whether it also 

leads to a loss of entanglement. In other words, could (R) also track the number of 

degrees of freedom of a system at scale R? which would imply (1.2) should be negative, i.e. 

(R) should be monotonically decreasing. 

For d = 2,^ a previous result of Casini and Huerta [6] shows that S2{R) is indeed mono- 
tonically decreasing for all Lorentz- invariant, unitary QFTs, which provides an alternative 
proof of Zamolodchikov's c-theorem [7]. 

In higher dimensions, the shape of S also matters. We argue in Sec. IV that5^P'^^'"^)(i?) 
has the best chance to be monotonic. At a fixed point, S^^^^'^'^\r) reduces to the previ- 
ously proposed central charge in all dimensions.'^ Its monotonicity would then establish the 
conjectured c-theorems [9-11] for each d. (For notational simplicity, from now on we will 
denote the corresponding quantities for a sphere simply as S{R) and Sd{R) without the 
superscript.) 

In Sec. V we consider a free massive scalar and Dirac field in ci = 3, where available 
partial results again support that^ 

SsiR) = Rq]^-S (1-3) 



^ for which S is given by two points and there is no need to have a superscript in 52 (i?). 

^ That the entanglement entropy could provide a unified definition of central charge for all dimensions was 

recognized early on in [8] and was made more specific in [10] including proof of a holographic c-theorem. 
^ A similar construction which involves the partition function (instead of entanglement entropy) has been 

used in [12] in connection with the g-theorem for d — 1. 
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is monotonic. 

In Sec. VI and VII we turn to holographic systems whose gravity dual satisfies the null- 
energy condition. In Sec. VI, among other things we show that when the central charges of 
the UV and IR fixed points are sufficiently close, Sd{R) is always monotonic in all dimensions. 
Sec. VII is devoted to numerical studies of various holographic systems in c? = 3 and (1 = 4. 
We find all the examples support the conjecture that: in d = 3, S^lR) is always non-negative 
and monotonically decreasing with R for Lorentz-invariant, unitary QFTs. 

In d = 4, where 



we find that while S^lR) appears to have the tendency to be monotonically decreasing, 
there exist holographic systems where it, however, is not always monotonic and can become 
negative. 

We conclude in Sec. VIII with a summary and a discussion of future directions. 
Note Added: When this paper is finalized, we became aware of [66] which has some overlap 
with our study. 

II. A REFINEMENT OF ENTANGLEMENT ENTROPY 

In our discussion below we will assume that the system under consideration is equipped 
with a bare short-distance cutoff 6q, which is much smaller than all other physical scales of 
the system. The continuum limit is obtained by taking 6o ^ while keeping other scales 
ffxed. The entanglement entropy for a spatial region is not a well-deffned observable in the 
continuum limit as it diverges in the 6q ^ limit. The common practice is to subtract 
the UV divergent part by hand, a procedure which is often ambiguous. The goal of this 
section is to introduce a refinement of the entanglement entropy which is not only UV finite, 
but also is most sensitive to the entanglement correlations at the scale of the size of the 
entangled region. 

A. Structure of divergences in entanglement entropy 

In this subsection we consider the structure of divergent terms in the entanglement en- 
tropy. We assume that the theory lives in fiat M^'"'^^ and is rotationally invariant. The 
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discussion below is motivated from that in [13] which considers the general structure of local 
contributions to entanglement entropy in a gapped phase. ^ We will mostly consider the 
vacuum state and will comment on the thermal (and finite chemical) state at the end. 

Let us denote the divergent part of the entanglement entropy for a region enclosed by a 
surface S as S'^^ . Then 5*^^^ should only depend on local physics at the cutoff scale near 
the entangling surface. For a smooth E, one then expects that S^j^ should be expressible in 
terms of local geometric invariants of S, i.e. 



where a denotes coordinates on S, F is a sum of all possible local geometric invariants 
formed from the induced metric hat and extrinsic curvature Kab of S. Note that here we are 
considering a surface embedded in flat space, all intrinsic curvatures and their derivatives 
can be expressed in terms of Kab and its tangential derivatives, thus all geometric invariants 
can be expressed in terms of the extrinsic curvature and its tangential derivatives. The 
proposal (2.1) is natural as S^J should not depend on the spacetime geometry away from 
the surface nor how we parametrize the surface. Thus when the geometry is smooth, the 
right hand side is the only thing one could get after integrating out the short-distance degrees 
of freedom. In particular, the normal derivatives of Kab cannot appear as they depend on 
how we extend S into a family of surfaces, so is not intrinsically defined for the surface itself. 

Here we are considering a pure spatial entangled region in a flat spacetime, for which 
the extrinsic curvature in the time direction is identically zero. Thus in (2.1) we only have 
Kab for the spatial normal direction. In more general situations, say if the region is not 
on a spatial hypersurface or in a more general spacetime, then E should be considered as 
a co-dimensional two surface in the full spacetime and in (2.1) we will have K^f^ with a 
running over two normal directions. 

Given (2.1), now an important point is that in the vacuum (or any pure state). 




(2.1) 



5'(s) ^ ^(2) 



(2.2) 



where S^^^ denotes the entanglement entropy for the region outside E, and in particular 



(2.3) 



div div ■ 



^ We thank Tarun Grover for discussions. 
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Recall that Kab is defined as the normal derivative of the induced metric and is odd under 
changing the orientation of S, i.e., in S^J it enters with an opposite sign. Thus (2.1) 
and (2.3) imply that F should be an even function of Kab- In a Lorentz invariant theory, 
there is also an alternative argument'' which does not use (2.2) or (2.3). Consider a more 
general situation with both K^j^ as mentioned above. The a index has to be contracted 
which implies that F must be even in K^^^. Then for a purely spatial surface we can just set 
the time component of K^j^ to zero, and F is still even for the remaining Kab- 

As a result one can show that for a smooth and scalable surface E of size R, the divergent 
terms can only contain the following dependence on R 

S!S = + «2i?"-' + • • • . (2.4) 

See Appendix A for a precise definition of scalable surfaces. Heuristically speaking, these 
are surfaces whose shape does not change with their size R, i.e. they are specified by a 
single dimensional parameter R plus possible other dimensionless parameters describing the 
shape. For such a surface, one can readily show that various quantities scale with R as (see 
Appendix A for more details) 

hab^R', Kab-^R, Da-^R^ (2.5) 

where Da denotes covariant derivative on the surface. As a result, any fully contracted 
quantity which is even in K, such as F in (2.1), can only give rise to terms proportional 
to with n a non-negative integer, which then leads to (2.4). Below we restrict our 

discussion to scalable surfaces. 

Now let us consider a scale invariant theory in the vacuum. On dimensional ground, the 
only other scale can appear in (2.4) is the short-distance cutoff ^o- We should then have 

11 , , 

«2~T^,--- (2.6) 

and so on. For odd d, the 0{R^) term is not among those in (2.4) and thus should be finite. 
For even d, there can be a log(5o term at the order 0(7?*^) and should come with log ^ in 
order to have to the right dimension. We thus conclude that for a scale invariant theory, 

^ We thank R. Myers for pointing this out to us. 
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the entanglement entropy across a scalable surface S in the vacuum should have the form 

+ . . . + A + e(^) + ^ + . . . odd d 

,.-2 + + So+^ ^) ' + R+ odd d 

|^ + --- + f + logf + const + | + --- even d 

where for notational simplicity we have suppressed the coefficients of non-universal terms. It 
is important to emphasize that S does not contain any divergent terms with negative powers 
of R in the limit 60 — )■ 0. The form (2.7) was first predicted from holographic calculations 
in [8] for CFTs with a gravity dual, s^^'* is an i?-independent constant which gives the 
universal part of the entanglement entropy. The sign factors before s)^ in (2.7) are chosen 
for later convenience. As indicated by the superscript, s)^ in general depends on the shape 
of the surface. 

For a general QFT, there could be other mass scales, which we will denote collectively as 
fi. Now the coefficients in (2.4) can also depend on /i, e.g., we can write ai as 

^1 = T^h^iM (2.8) 

and similarly for other coefficients. Note that by definition of 6q, we always have /i^o ^ 1 
and hi can be expanded in a power series of /i(5o. Now for a renormalizable theory, the 
dependence on /i must come with a non-negative power, as when taking fi6Q — )■ 0, ai should 
not be singular and should recover the behavior of the UV fixed point. In other words, for 
a renormalizable theory, the scale (s) /i arises from some relevant operator at the UV fixed 
point, which implies that fiSo should always come with a non-negative power in the limit 
^6q — > 0. This implies that the UV divergences of Oj should be no worse than those in (2.6). 
In particular, there cannot be divergent terms with negative powers of R for even d, and for 
odd d the divergence should stop at order 0{R). These expectations will be confirmed by 
our study of holographic systems in Sec. VIC and VI D (see e.g. (6.26)), where we will find 
that /ii(/i5o) has the expansion hi{fi6o) = cq + C2(/i5o)^" + C3(/i5o)^" + ■ ■ ■ , where a = d — A 
with A the UV dimension of the leading relevant perturbation at the UV fixed point. 

So far we have been considering the vacuum. The discussion of the structure of diver- 
gences should work also for systems at finite temperature or finite chemical potential. In 
such a mixed state, while (2.2) no longer holds, equation (2.3) should still apply as the short- 
distance physics should be insensitive to the presence of temperature or chemical potential. 
Also recall that for a Lorentz invariant system, there is an alternative argument for (2.4) 
which does not use (2.3). 



B. Properties of 



Given the structure of divergent terms in S^^\R) discussed in the previous subsection, 
one can then readily check that when acting on S^^\R) with the differential operator in (1.1), 
all the UV divergent terms disappear and the resulting (R) is finite in the continuum 
limit 6q ^ 0. In fact, what the differential operator does is to eliminate any term (including 
finite ones) in S^^^R) which has the same /^-dependence as the terms in (2.4). We believe, 
for the purpose of extracting long range correlations, it is sensible to also eliminate possible 
finite terms with the same /^-dependence, as they are "contaminated" by short- distance 
correlations. In particular, in the continuum limit this makes iS^^''(i?) invariant under any 
redefinitions of the UV cutoff 6o which do not involve R.^ With a finite 6o, (R) does 
depend on 6o, but only very weakly, through inverses powers of ^. This will be important 
in our discussion below. 

In the rest of this section we show that the resulting (R) is not only UV finite, but 
also have various desirable features. In this subsection we discuss its behavior in the vacuum, 
while in Sec. II C discuss its properties at a finite temperature and chemical potential. 

For a scale invariant theory, from (2.7) we find that for all d 

Sf\R) = sP (2.9) 

is i?- in dependent. The sign factors in (2.7) were chosen so that there is no sign factor 
in (2.9). Note that if we make a redefinition of the form 6o — > 5o (1 + cifi6o + C2(/i(5o)^ + ■ ■ ■) 
where fi is some mass scale, for odd d the UV finite term in (2.7) is modified. But (-R), 
and s^^'' as defined from (2.9), is independent of this redefinition. 

Let us now look at properties of SP{R) for a general renormalizable QFT (i.e. with a 
well-defined UV fixed point). Below we will find it convenient to introduce a fioating cutoff 
6, which we can adjust depending on scales of interests. At the new cutoff 6, the system is 
described by the Wilsonian effective action leffi^', ^o), which is obtained by integrating out 
degrees of freedom from the bare cutoff Sq to S. The entanglement entropy S^^\R;So,S) 
calculated from leffi^', ^o) with cutoff 6 should be independent of choice of 6. So should the 
resulting Sj^\R). Below we will consider the continuum limit, i.e. with bare cutoff 6^ — > 0. 

® Since R is the scale at which we probe the system, reparameterizations of the short-distance cutoff should 
not involve R. 
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First consider the small R limit, i.e. R is much smaller than any other length scale of 
the system. Clearly as i? — > 0, these other scales should not affect S^'^^R), which should be 
given by its expression at the UV fixed point. Accordingly, (R) also reduces to that of 
the UV fixed point, i.e. 

5f ^(i?) ^ R^O. (2.10) 

As we will see in Sec. VI E, studies of holographic systems (with S given by a sphere) predict 
that the leading small R correction to (2.10) is given by 

S^P{R) = '^^^ + 0((/ii?)2°) R-^0 (2.11) 

where a = d — A with A < d the UV dimension of the leading relevant scalar perturbation. 
Equation (2.11) has a simple interpretation that the leading contribution from a relevant 
operator comes at two-point level. We believe it can be derived in general, but will not 
pursue it here. 

The story is more tricky in the large R limit, as all degrees of freedom at scales between 
the bare UV cutoff 6o and R could contribute to the entanglement entropy S^'^^R) in this 
regime. Nevertheless, one can argue that 

SPiR) ^ sf'^^'l R^oc (2.12) 

as follows. When R becomes much larger than all other length scales of the system, we can 
choose a floating cutoff 6 to be also much larger than all length scales of the system while 
still much smaller than R, i.e. 

— (5 (2.13) 

Ail fJ'2 

where fii,i = 1,2, ■■■ denote possible mass parameters of the system. Now the physics 
between 6 and R is controlled by the IR fixed point, i.e. we should be able to write 
S^^\R) again as (2.7), but with 6o replaced by 6, and s^^^ by s^J^'^^K Then equation (2.12) 
immediately follows. In other words, while in terms of the bare cutoff 6o, the entangle- 
ment entropy S^^\R) = S'-^\6o,R, /^i;/^2,''') could be very complicated in the large R 
regime, involving many different scales, there must exist a redefinition of short-distance cut- 
off = 5ir((5, /ii, /i2, ■ ■ ■ ), in terms of which S^^\R) reduces to the standard form (2.7) 
with 6 replaced by (5ir, and s^^"* by s^^'^^\ In fact, higher order terms in (2.7) with negative 
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powers of R also imply that generically we should expect the leading large R corrections 
to (2.12) to have the form 

5f)(fl) = ,f fl^oo. (2.14) 

[O(^) evend 

This expectation is supported by theories of free massive scalar and Dirac fields as we will 
see in Sec. V, and by holographic systems as we will see in Sec. VI E. Holographic systems 
also predict an exception to (2.14) which happens when the flow away from the IR fixed 
point toward UV is generated by an irrelevant operator with IR dimension Air sufficiently 
close to d, for which we have instead (see Sec. VI E)^ 



a < ^ odd d 
a < 1 even d 



'^f (^) = 4'''^^ + 0(^^), for <! ' , i?^oo, (2.15) 



where a = Air — d. 

By adjusting the floating cutoff 6, one can also argue that (i?) should be most sensitive 
to contributions from degrees of freedom around R. Consider e.g. a length scale Li which is 
much smaller than R. In computing S^^\R), we can choose a floating short-distance cutoff 
6 which satisfies 

Li<^5<^R. (2.16) 

As discussed at the beginning of this subsection, by design S^\R) is insensitive to short- 
distance cutoff 6 when 6 <ti R.^^ We thus conclude that sI^\r) should be insensitive to 
contributions of from d.o.f around Li. 

While our above discussion around and after (2.12) assumes a conformal IR fixed point, 
the discussion also applies to when the IR fixed point is a gapped phase, where there are 
some differences depending on the spacetime dimension. For odd d, using = 3 as an 
illustration, the entanglement entropy for a smooth surface S in a gapped phase has the 



^ The expression below is derived in Sec. VIE for E given by a sphere and closely separated UV/IR fixed 
points. We believe the result should be more general, applicable to generic systems and smooth E, but 
will not pursue a general proof here. 
^'^ Of course ultimately as mentioned earlier sj^\R) should be independent of choice (5, when one includes 
all possible dependence on S including those in coupling constants. Here we are emphasizing that even 
explicit dependence on S should be suppressed by negative powers of 
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form (see e.g. also [13]) 

S^^\R) = aR - + 0{R-^) (2.17) 
where 7 is the topological entanglement entropy [14, 15]. We then have 

SP{R)^-f, R^oo. (2.18) 

In gapped phases without topological order, 7 = 0. Thus a nonzero s!^\r — )■ 00) signals the 
system has long range entanglement, i.e. the system is either gapless or topological-ordered 
in the IR. The two cases can be distinguished in that for a topological ordered phase 7 
should be shape- independent, but in a gapless case, s'"^'^^^ in (2.10) is shape-dependent. 

For even d, in a gapped phase we expect that S^'^\R) does not have a term proportional 
to log R for large R, and thus we should have 

s£\R)^0, R^oo, n = l,2,---. (2.19) 

Nevertheless, it has been argued in [13] that the size-independent part of the entanglement 
entropy contains topological entanglement entropy. Such a topological term could not be 
captured by (R), as all terms in (1.1) contain derivatives with respect to R for even d. 
This is not surprising, as in even d, the i?-independent part of the entanglement entropy also 
contains a finite non-universal local part, as is clear from the discussion around (2.4). Thus it 
is not possible to separate the topological from the non-universal contribution using a single 
connected entangling surface, and one has to resort to constructions like those in [14, 15] to 
consider combination of certain regions in such a way that the local part cancels while the 
topological part remains [13]. 

C. Finite temperature and chemical potential 

As discussed at the end of Sec. II A, we expect (R) should also be UV finite in the 
continuum limit at a finite temperature or chemical potential. Here we briefly discuss its 
properties, and for simplicity will restrict to a scale invariant theory. 

See also [16] regarding scaling behavior of the entanglement entropy at finite T. In ref. [16] considered 
the entanglement entropy itself with UV part subtracted manually. 
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For a scale invariant system at a finite temperature T, since there is no other scale in the 
system, Sj^\R,T) must have a scaling form, i.e. 

S^f\R,T) = Sf\RT) . (2.20) 

In particular, in the high temperature limit, i.e. RT ^ 1 it must be dominated by thermal 
entropy at leading order, while in the low temperature limit RT — )■ 0, it should reduce to s)^^ . 
For a scale invariant theory, the thermal entropy has the form Sj^ = rjaT'^'^Vj] ~ (TRY~^, 
where Vs is volume of the spatial region enclosed by S and rjd is some constant. Thus we 
should have 

rv^ f-^r oddd 

SPiRT)^^ , RT^oo. (2.21) 

[[fc§T'^r evend 
More explicitly, for d = 3 we expect when RT ^ 1, 

(RT) = r/,TVs + cf ^ + ■ ■ ■ (2.22) 

where the first term is simply the thermal entropy, ■ ■ ■ denotes terms with negative powers 
of RT. The second term q Ms a constant. It would be interesting to compute this constant 
for some explicitly examples to see whether some physical interpretation (or significance) 
can be attached to it. 

Similarly, with a nonzero chemical potential generalization of (2.20) we expect 

that 

5f (i?,/i,r)=5f)(^^,i?r^ . (2.23) 

III. ENTANGLEMENT ENTROPY OF A (NON)-FERMI LIQUID 

In this section we show that the entanglement entropy of a (non)-Fermi liquids can be 
obtained by simple dimensional analysis. Consider a rf- dimensional system of a finite 
fermions density whose ground state is described by a Fermi surface of radius kp. We 
have in mind a Fermi liquid, or a non-Fermi liquid described by the Fermi surface coupled 
to some gapless bosons (as e.g. in [21]). In either case, the low energy dynamics of the 



See also a recent discussion in [16] based on finite temperature scaling and crossover. See also [17-20] for 
recent discussion of logarithmic enhancement of holographic "non-Fermi liquids." 
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system involves fermionic excitations locally in momentum space near the Fermi surface, 
and different patches of the Fermi surface whose velocities are not parallel or anti-parallel 
to each other essentially decouple. In particular, kp drops out of the low energy effective 
action. We thus expect in the large R limit, the "renormalized entanglement entropy" 
S^\r) should be proportional to the area of the Fermi surface Aps oc which can be 
considered as the "volume" of the available phase space. Since there is no other scale in the 
system than R, (R) should then have the form 

S^iR) oc kp^R^'^ oc ApsA^, R-^oo . (3.1) 

where A^, denotes the area of the entangling surface S. In other words, our "renormalized 
entanglement entropy" should satisfy law." Using (1.1) one can readily see that 

the area law (3.1) translates into the well-known behavior in the original entanglement 
entropy [22, 23] (see also [24]) 

S^{R) oc kp^R''-^ \og{kpR) + ■ ■ ■ oc ApsA^ log(A^s^s) + • • • , (3.2) 

where kp in the logarithm is added on dimensional ground and ■ ■ ■ denotes other non- 
universal parts. We note that this result does not depend on whether the Fermi surface has 
quasi-particles or not, i.e. whether it is a Fermi or non-Fermi liquid, only depends on the 
expectation that iS^ (i?) is proportional to the area of the Fermi surface. 

This analysis can also be immediately generalized to predict the qualitative behavior of 
the entanglement entropy of higher co- dimensional Fermi surfaces. For a co-dimensional n 
Fermi surface we should have^^ 

sf\R) oc {kpRf-'' (3.3) 

which implies that in the entanglement entropy itself^^ 

f ikpRy-'' \ogikpR) n even 
S^{R) oc r . (3.4) 

UkpRy-"" nodd 

Thus we find that there is a log R factor only for even co-dimensional Fermi surfaces. These 
results are again independent of whether there are quasi-particles. Note that for a Fermi 



We define the co-dimension with respect to the fuU spacetime dimension d. 
These results were also obtained by B. Swingle (unpublished). 
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point where n = d, equation (3.4) in consistent with one's expectation that for massless 
fermions there is a universal logR term only for even d. For general n, at least for free 
fermions, the alternating behavior of logarithmic enhancement in (3.4) may also be under- 
stood (by generalizing an argument of [24]) as follows: at each point of a co-dimensional n 
Fermi surface, there is an n-dimensional free fermion CFT. The logi? appearance in (3.4) 
is then consistent with the fact that for an n-dimensional CFT, there is a universal logi? 
piece only for n even. 

It would be interesting to see whether our discussion may also be used to understand 
the logarithmic enhancement in the entanglement entropy of the critical spin liquids in [25] 
which are described by a projected Fermi sea state. 

IV. ENTANGLEMENT ENTROPY AS MEASURE OF NUMBER OF DEGREES 
OF FREEDOM 

For the rest of this paper we will restrict our discussion to the vacuum. In Sec. II we 
showed that in the vacuum (i?) introduced in (1.1) has various desirable features: 

1. It has a well-defined continuum limit. 

2. For a CFT, it is independent of R and given by the universal part of the entanglement 
entropy (2.7) 



3. For a renormalizable quantum field theory, it interpolates between the values of UV 
and IR fixed points as R is increased from zero to infinity. 



It should be understood that in (4.2) if the IR fixed point is described by a gapped 
phase, then is either given by the topological entanglement entropy (for odd d) 
or zero (for even d). 

4. It is most sensitive to degrees of freedom at scale R. 



(4.1) 




(4.2) 



16 



Thus (R) provides a set of observables which can be used to directly probe and char- 
acterize quantum entanglement at a given scale R. As discussed in the Introduction, these 
properties also imply that we may interpret the dependence on i? as a RG flow. A natural 
question which then arises is whether iS^^'' (R) could also provide a scale-dependent measure 
of the number of degrees of freedom of a system. Given the physical intuition that RG leads 
to a loss of degrees of freedom, a necessary condition for this interpretation is then 



which in turn requires (given (4.2)) 



^dSf\R)^ < (4.3) 
dR 



4"'""' > 4=^'"' . (4.4) 



Note that (4.4) alone is enough to establish s)^ as a measure of the number of degrees of 
freedom for CFTs, while establishing (R) as a measure of degrees of freedom for general 
QFTs requires a much stronger condition (4.3). 

For d = 2, the entangled region becomes an interval (there is no shape difference) and 
equation (1.1) reduces to^^ 



which for a CFT then gives [26-28] 



S2{R) = R^ , (4.5) 



^2 = ^ (4.6) 



where c is the central charge. In this case, Zamolodchikov's c-theorem [7] ensures (4.4) 
and there exists a beautiful proof by Casini and Huerta [6] showing that S2{R) is indeed 
monotonically decreasing for all Lorentz- invariant, unitary QFTs. Note that while there 
already exist an infinite number of c-functions [30] including Zamolodchikov's original one, 
S2{R) has some special appeal, given that it also characterizes the entanglement of a system. 
We would like to propose that it gives a "preferred" c-function which best characterizes the 
number of d.o.f. of a system at scale R. 

In higher dimensions, the shape of S also matters. Could (4.3) and (4.4) apply to generic 
or only certain shapes? For this purpose, consider first the weaker condition (4.4). 

The function (4.5) has also been discussed e.g. in [5, 29] as the universal part of the entanglement entropy 
in d = 2. 
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For even d > 4, since s\ appears as the coefficient of the divergent term log^o in (2.7), 
it can be expressed in terms of integrals of the geometric invariants associated with S (re- 
call (2.1)), and in particular related to trace anomaly [8, 31]. For example, for o? = 4 [8, 31] 

sf) = 2a4 [ d^aVhE2 + CA I d^aVhh (4.7) 
where 04 and C4 are coefficients of the trace anomaly^^ and (below K = Kabh"'^) 

E2 = ^n=^ [K^ - KabK'^^) (4.9) 
47r 47r 

/2 = ^ {^-K^ - KabK'^'^^ . (4.10) 

In (4.9) TZ is the intrinsic curvature on S and in the second equality of (4.9) we have used 
the Gauss-Codacci relation in flat space. E2 is the Euler density for d = 2 and d'^a\/hE2 
is a topological invariant with value 2 for a surface with spherical topology. I2 is a Weyl 
invariant and is zero for a sphere. For a sphere we then have 

^(sphere) ^ (4.11) 

while for other shapes, will be a linear combination of and C4. More than twenty 
years ago, Cardy conjectured that [9] 04 and its higher dimensional generalizations obey the 
analogue of c-theorem. Only very recently was it proven for d = 4 [32]. In addition, there 
are strong indications any combination of 04 and C4 (including C4) will not satisfy such a 
condition [32]. Thus for d = 4 only for E = sphere could the condition (4.4) be satisfied. 
For higher even dimensions the situation is less clear, but one again has [8, 10, 33]^^ 

^(sphere) ^ _ ^^^^2) 

For odd d, s\ does not arise from local terms in (2.7), thus we do not expect that it 
can be expressed in terms of local geometric invariants on S. It was found in [33] that for 



We use the convention 

{T\) = -^WabcdW^'"' + 2aiEi (4.8) 

with Wabcd the Weyl tensor and = ^ {RabcdR"-^'"^ ~ -iRabR"^ + R^) the Euler density. The relation 
to A and B of [31] is ^ = B = Also note there is a minus sign in the definition of s"l^^ in (2.7). 
In comparing with [31] we have also set the extrinsic curvature in the time direction to zero. 
In addition, the structure of (4.7) persists in higher even dimensions except that there are more Weyl 
invariants [10, 34]. 
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a sphere 

4^'''^^^^ = (l0g^)finite (4.13) 

where (logZ)finite is the finite part of the Euchdean partition for the CFT on 5"^. 
Some support has been found that s^'^p'^'^''*^^ (equivalently (log Z)finite) satisfies the condi- 
tion (4.4) [10, 11]. 

To summarize, for E given by a sphere, there are (strong) indications that s^'^p*^'^''^^ could 
satisfy (4.4) and thus provide a measure of the number of degrees of freedom for CFTs in 
all dimension (including both odd and even) [10, 33]. Below we will simply call s^^'^p'^'^''*'^ 
the central charge of a CFT. For (i 7^ 4, other shapes could still provide a similar measure, 
which will be left for future investigation. For the rest of the paper we study the stronger 
condition (4.3) for E = sphere in c? = 3 and d = 4. For notational simplicity, we will drop 
the superscript "sphere" in various quantities and denote them simply as S{R),Sd{R) and 

Sd- 

V. FREE MASSIVE SCALAR AND DIRAC FERMIONS IN d = 3 

In this section we consider S^lR) for a free massive scalar and Dirac field. For a free 
massive field, we expect that S^lR) should approach that for a massless field as i? — )■ 0, 
and as — 7- 00. We would like to see whether it is monotonic and positive in between. 
Recently, in the limit of mR 3> 1, it was found in [5, 35] that {S is a short-distance cutoff) 

5sealar(mi?) = # f " ^ " ^ + " " " (5-1) 

5Dirac(mi?) = #:^-^mi?-^^ + -- -. (5.2) 
40 mR 

From (1.3) we thus find that 

Wmfl) = + (5.3) 

So„„(mfi) = +|,^ + --- (5.4) 

which are indeed monotonically decreasing with R. Note that the ^ fall-off in the above 
expressions is also consistent with our earlier expectation (2.14). 

We emphasize that if one simply subtracts the divergent part in (5.1), then the resulting 
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approaches minus infinity linearly as i? — )■ oo and thus does not have a good asymptotic 
limit. The presence of such a linear term can be understood as a finite renormalization 
between the short distance cutoffs of the UV and IR fixed points, as discussed in Sec. II. In 
contrast, Ssi^R) approaches zero as -R — )■ oo as one would expect of a system with a mass 
gap. 

We have also calculated S^lR) numerically for a massive scalar field^^ for a range of mR 
as shown in Fig. 1 (see Appendix B for details for the numerical calculation). The numerical 
result is consistent with our expectation of the limiting values of S^lR) in the small and 
large R limits, and also suggests that it is monotonic in between. 

S3(K) 




J , , , , _L , , , , _L , , , , _L , , , , _L , , , , _L , , , , 

FIG. 1: 53(i?) for a free massive scalar field: The red point is the value for m = 0. The black dashed 
line is the result of the asymptotic expansion (5.3). The numerical results are computed following 
the method of [4] with a radial lattice discretization. We choose the system size to be Lj/j = 200a, 
where a is the lattice spacing. To avoid boundary effects the restriction to 10a < R < 45a was 
made. To extend the range of mR we obtained the results for 1/m = 20a, 40a, 120a. In the plots, 
the orange dots are data points for 1/m = 120a, the blue ones are for 1/m = 40a, and the green 
ones are for 1/m = 20a. As expected all our data points collapse into one curve as 53(i?) can only 
depend on mR in the continuum limit. For more details see Appendix B. 



Compared to that of a scalar, the computation for a Dirac fermion requires significantly more computer 
time to achieve the same accuracy. We will leave it for future investigation. 
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VI. Sd{R) FOR HOLOGRAPHIC FLOW SYSTEMS 



In this section we discuss properties of Sd^R) (defined for a sphere) for systems with a 
gravity dual using the proposal of [8, 36] (see [37] for a review), which relates the entan- 
glement entropy to the area of a minimal surface. Other recent discussion of entanglement 
entropy in holographic RG flow systems include [10, 38-42] 

We will restrict our discussion to > 3. After a brief discussion of the general set-up, we 
derive a relation between Sd and the undetermined constant in the asymptotic expansion of 
the minimal surface. We then show that when the central charges of the IR and UV fixed 
points are close, for all dimensions, is always monotonically deceasing with R at leading 
order in the expansion ot difference of central charges of the UV and IR fixed points. Thus 
for flows between two sufficiently closely separated fixed points, S{R) appears to provide a 
good central function. 



A. Gravity set-up 

We consider a bulk spacetime which describes a renormalization group flow in the bound- 
ary theory. We assume that the system is Lorentz invariant. The flow can be induced either 
by turning on the source for some relevant scalar operators or by certain scalar operators 
developing vacuum expectation values (without a source). Below we will refer to them as 
source and vev deformation respectively. We denote the corresponding bulk fields by (pa- 

The bulk action can be written as 



1 

2^ 



R - lGabd(j)ad(j)b - V{<f)a) 



(6.1) 



where Gab is some positive-definite metric on the space of scalar fields. The spacetime metric 
can be written in a form 

ds^ = — i-dt^ + dp^ + p^dnl_, + — \ . (6.2) 

We assume that V{(j)a) has a critical point at 0a = with V{0) = — ^^^r^, which corresponds 
to the UV fixed point. Near the boundary z = 0, 

0a(z)->O, /(^)^1, ^^0 (6.3) 
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and the spacetime metric is that of AdS^+i with curvature radius L. Einstein equations and 
positive-definiteness of the kinetic term coefficients Gab require that the evolution of / with 
z should satisfy [43] 

dj{z) > (6.4) 

i.e. f{z) is a monotonically increasing function. More generally, equation (6.4) follows from 
the null energy condition and Einstein equations, regardless of specific form of the scalar 
action. 

At small z, f{z) can be expanded as 

f{z) = 1 + /i2°;z2a ^ . . . ^ ^ ^ q ^g g^ 

where /i is some mass scale and a some positive constant. For a source deformation, a = 
d — A with A the UV dimension of the leading relevant perturbating operator (i.e. the one 
with the smallest a).^^ For a vev deformation we have a = A."^^ 
As z — )■ oo, we can have the following two possibilities: 

1. Flow to an IR conformal fixed point. In this case (pa approaches a neighboring critical 
point (f)a = (t)T with ^(0^°)) = and \/(0(°)) < 1/(0), L > Ljr. The fiow 

^IR 

solution then describes a domain wall with the metric (6.2) interpolating between two 
AdSd+i with curvature radius given by L and Ljr respectively, i.e. 

f{z) > 1, z ^ +00 . (6.6) 

^IR 

Near the IR fixed point, i.e. z — > oo, / can be expanded as 

1_ 

l]r if^^y 

where a = Air— c/ with Air the dimension of the leading irrelevant perturbing operator 
at the IR fixed point and /i is a mass scale characterizing irrelevant perturbations. 

2. Flow to a "gapped" phase: since / has to increase with z, instead of approaching a 
constant as in (6.6), f{z) can blow up as z — )■ oo, e.g. 

f{z) oc 2", n>0 . (6.8) 



f(') = T^-77-^ + --- (6-7) 



Note that when A = | for which a = |, we should replace the second term in (6.5) by {fizY {log ^z)"^ . 
The above description is for the standard quantization. In the alternative quantization (which applies to 
| — 1<A< |), we have instead a = A for a source deformation and a = d — A for a vev deformation. 
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The spacetime then becomes singular at z = oo. Given n > 0, the z = oc singularity 
in fact sits at a finite proper distance away. From the standard IR/UV connection 
we expect that the system should be described by a gapped phase in the JR. Exphcit 
examples include the GPPZ [44] and Coulomb branch flow [45, 46] which we will 
discuss in more detail in next section. While one should be wary of such singular 
spacetimes, they appear to give sensible answers for correlation functions consistent 
with the interpretation of a gapped phase (see e.g. [47]).^^ In this paper we will assume 
such singular geometries make sense. 



B. Holographic entanglement entropy 

From the prescription of [36], the entanglement entropy for a spherical region of radius 
R is obtained by 

O-rrTd—l n Td~l 

S{R) = ^uj,^2A = KA, K= — a;,_2 (6.9) 

where uJd-2 is the area of a unit {d — 2)-dimensional sphere and A is obtained by minimizing 
the surface area 

rz^ d-2 j ]~ 

with the boundary condition at infinity 

p[z = Qi) = R. (6.11) 

Depending on the spacetime metric, there can be two kinds of minimal surfaces as indicated 
in Fig. 2. For the disk type, the minimal surface ends at a finite Zm with 

p{Zm) = 0, p{Zm) = OO . (6.12) 

The cylinder type solution extends all the way to 2; = oo with 

p{z) -)■ po, p'iz) -)■ 0, z 00 (6.13) 



In Coulomb branch flow [45, 46] which describes a Higgs phase of the Af = 4 SYM theory, there is a single 
Goldstone mode corresponding to spontaneous breaking of conformal symmetry. In the large N limit, the 
effect of such a gapleess mode on observables like entanglement entropy can be neglected. We will thus 
still call it a "gapped" phase. 
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with po a finite constant. When the IR geometry is given by another AdS^+i as in (6.6), 
the minimal surface is always of disk type. In Appendix C we show that for singular IR 
geometries with (6.8), cylinder type solution can exist when the exponent n > 2. As also 
discussed in Appendix C, for n > 2, there exists a critical minimal surface which closes off 
exactly at 2; = 00, i.e. it can be viewed either as a surface of cylinder topology with po = or 
a surface of disk topology with Zm = 00. This critical surface will be important for "second 
order phase transitions" in entanglement entropy we observe in some examples of Sec. VII. 



2 — > 00 




* p 



FIG. 2: Cartoon of a minimal smface of disk topology (black) v.s. a minimal surface of cylinder 
topology (red). The cylinder type surface is possible only for (6.8) with n > 2. 

At the UV fixed point, / = 1, the minimal surface is given by the hemisphere [8] 

P0(Z) = - Z^, ZmO = R. (6.14) 

Evaluating on (6.14), the action (6.10) is divergent near z = 0. We introduce a UV cutoff 
z = 6 and the regularized action becomes 

■4 = ^/"<i^^^^^3#^- (615) 

One then finds the entanglement entropy has the form (2.7), with 

<,(uv) _ (UV) _ (^-3)!! , . 
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Similarly at the IR fixed point with / = -M- we find 



o(iR) (IR) 27rLjj^^ /C1^7^ 

^ (rf-2)!! ' ^ ^2 ^^^-2 • (6-17) 

C. 5d in terms of asymptotic data 

We now derive a formula to relate Sd directly to a finite coefficient in the asymptotic 
expansion of the minimal surface solution near the boundary. For definiteness we will restrict 
our discussion to a single scalar field, but the conclusion is general. 

Consider the action (6.10) with a cutoff ai z = 6, 



A = I ^ dz^Jp'^ + \ = I dzC (6.18) 



15 ^ V J J5 

with boundary conditions (6.11) and (6.12) (or (6.13) with Zm = oo). In the above equation, 
varying R with 6 fixed, and using the standard Hamilton- Jacobi method we find 



dR~ ^^"^'dR dR 



Ii = ^ = P^^ / , V, = Iip'-C = - P ^ / : (6.20) 



where 



/ J ' f 

and we have used that 'H{zm) = due to the boundary condition there (for both (6.12) 

and (6.13)). Equation (6.19) only depends on the local solution near infinity and can be 

evaluated using the asymptotic expansion there. 

For small z, f{z) has the expansion (see e.g. also [48]) 

oo oo 

f{z) = 1 + 5^ /n^'+" + Yl C"^"" (6-21) 

n=0 m=2 

where a was introduced in (6.5). Note that here we are considering a fiat boundary, so the 
first series of (6.21) starts only at order 0{z'^). The second series comes from the scalar 
contribution. Similarly, p{z) can be expanded as [48-50]^^ 

2R 



p{z) = R-^ + --- + UR)z'' + ---+ J2 «n™(i?)^"+"^" . (6.22) 



n=2,m,=2 



There are also logarithmic terms when two series share common terms. For example, when n + ma = d, 
there is a term proportional to z"^ log fiz with fi a mass scale associated with relevant perturbation. 
Presence of these additional logs will not affect our discussion. So below we will not make them explicit. 
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In (6.22), the first ■ ■ ■ contains only even powers of z up to z'^ and tlie second ■ ■ ■ contains 
integer powers greater tlian d. Note tliat tfie structure of expansion in (6.22) excluding the 
last term follows from the standard Fefferman-Graham expansion applied to a submanifold in 
AdS [49], while the last term comes from the similar series in the expansion of /. In (6.22), q 
is undetermined and all other coefficients can be determined in terms of coefficients in (6.21) 
and Cd from local analysis of the minimal surface equation of motion near 2; = 0. 
Plugging the expansion (6.22) into (6.19) we find that in the limit 5—7-0 



1 - -''^-'-^^^ 



(6.23) 

non— universal 



where 

.4(l.(-in(l-(-l)^)|f^ (6,24) 

in non-vanishing only in = 4A;, k = 1, 2, ■ ■ ■ dimensions. The non-universal part has the 
form 

^ A 

= bi{5)R'^-^ + b2i5)R'^-^ + --- (6.25) 



dA 
dR 



non— universal 

with only terms of non-negative powers of R. bn can be expressed in terms of a divergent 
series of 5, with 5'^~'^ the most divergent term in 61, 5^"'^ for 62? and so on. More explicitly, 
introducing a dimensional parameter /i to characterize the dimension of c„ (with c„ ~ /i""), 
then the series in bi can be schematically written as 

61 ~ ^ (1 + (M)'" + (/i5)'" + ■■■)+ const . (6.26) 

The structure of non-universal terms (6.25) are consistent with that argued in Sec. II (recall 
the discussion below (2.8)). They go away when acted with the differential operator in (1.1) 
and do not appear in Sd- 

As a result Sd{R) can be solely expressed in terms of Cd{R)- More explicitly, from (6.23) 
and (1.3) we find that for d = 3 

^So(R^ = -3n^r,jm + 3 / 

where C is determined by requiring iS3(-R = 0) reduces to the value at the UV fixed point. 
Also note 

^ "^^^ -SdniRcsm . (6.28) 



1 

—SsiR) = -3R^C3{R) + 3 / dRRcsiR) + C (6.27) 
K Jo 



K dR 
Similarly, for = 4, 



= 1 - 2R:'ci{R) - 2R^^ . (6.29) 
K dR 
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D. Two closely separated fixed points 



Let us now consider the situation in which the bulk cosmological constants of the UV 
and IR fixed points are close to each other, which translates into the boundary theory as 
that the central charges Sd (6.16) of the UV and IR fixed points are close. In this case we 
have 




f{z) = l + eg{z), e= g'iz)>0, giz) ^ { . (6.30) 

-^IR 

We expect the minimal surface to be close to that of a CFT, i.e. 

p(z) = po{z) + epi{z) + ■ ■ ■ (6.31) 

with po given by (6.14), and 

A = Ao + eAi + -- - . (6.32) 
Under variations of /, the variation of the cutoff action (6.18) can be written as 

5A = dz—6f-nizm)5z^-mpl (6.33) 

where H and 11 were introduced in (6.20) and Szm and Sp denote the induced variations due 
to 6f. Now apply the above equation to (6.30)-(6.32), we find that 

A= dz— giz) + ni5)\^^^p,i5) (6.34) 

where T-L{zm) again vanishes due to the boundary condition (6.11) at Zm- Pi{z) is obtained 
by solving the linearized equation following from (6.10) around (6.14). To find its value at 
the short-distance cutoff 6, it is enough to use the leading order expression (6.5) for / (with 
g{z) = (pz)"^" + ■■■)' we find that^^ 

+ e(a) ^ ^(^^i±|^|-^ . (6.35) 

Plugging the explicit form (6.14) of po ^^^d (6.35) into (6.34), we find that 

IS 



Ai = -^£dz (R' - z'Y-^ z'-'g{z) + c{a) ' (/i^)^" + ■ ■ ■ • (6.36) 



23 When a = equation (6.35) becomes pi{S) = -^(M'^)''"^ yjprp logAt-^. 
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The divergent terms in the above equation are again consistent with the structure described 
in Sec. II A and thus will not contribute to ^^(i?).^^ 
Now applying (1.1) to (6.36) we find that, for odd d, 

-(l),o^_ (<i-l)!!A' (d-iy.\K " 







and for even d 



r^^^-^p^'-^ r^^^. (6.38) 



{d-2)\\2RJo VR^ - {d-2)U2Jo yT^x 

It is interesting to note that for a monotonic g{z), both (6.37) and (6.38) are monotonic 

^ < . (6.39) 

Also note that in the limit R oo, we can replace g{xR) in both expressions by its value 
at infinity g{oo) = 1, leading to 

which is precisely the expected difference between the values of the IR and UV fixed point 
(see (6.16) and (6.17)) expanded to first order in e. Thus to first order in e we have recovered 

S,{R -> oo) = . (6.41) 
E. Leading small and large i?-dependence 

In this section we examine the leading small and large i?-dependence of Sd{R)- 

When R is small, the range of z covered by the minimal surface is also small. Thus we can 

use the asymptotic form (6.5) of the function f{z) = 1 + 5/ H with 6f = (/iz)^", and treat 

the deviation from pure AdS as a small perturbation. To first order in 6f, equation (6.34) 
then again applies with g there replaced by 6f. The calculation then becomes very similar 
to that of the last subsection and we find the following leading behavior 

5,(i?) = 5r^-|^^|a(a)(/ii?)2" + ---, R^O (6.42) 

Note that for some special values of a, both terms in the above equation can lead to logarithmic terms 
which contain log fid ^ as pointed out before in [5, 48]. These terms are also consistent with the structure 
described in Sec. II. 
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with 5^^^^ given by (6.16) and a{a) is given by 



1 



dodd 



a(a) = { . (6.43) 

deven 

Recall that for a source deformation a = d — A where A is the UV dimension of the 
leading relevant operator. Equation (6.42) has a simple interpretation that the leading 
contribution from a relevant operator comes at two-point level, and should be valid also 
outside holographic systems. Also note that the sign of the numerical coefficient in (6.42) is 
such that at this order Sd{R) is monotonically decreasing with R. 

In the large R limit, estimating the leading i?-dependent correction to the asymptotic 
value Sd{R — > oo) = Sf^^ for a general holographic system becomes more challenging. Here 
we will consider the case of two closely separated fixed points discussed in the last subsection 
where we can take advantage of the closed expressions (6.37) and (6.38). At large R, the 
argument xR for g{xR) is large except near x = 0, thus we should be able to use the 
asymptotic expansion of g{z) for large z (from (6.7) and (6.30)), 

g{z) = 1 - ifiz)-^^ + ■ ■ ■ . (6.44) 

Recall that a = Air — d where Air is the IR scaling dimension of the leading irrelevant 
operator at the IR fixed point. Plugging (6.44) into (6.37) and (6.38), the 1 in (6.44) will 
give (6.40). To to first order in e we can write 

5,(/2)=5r^ + J,(i?) + 0(a (6.45) 

where 

((] _ eK 1 IH odd d 

In (6.46) the ■ ■ ■ denotes contributions from higher order terms in (6.44) (represented by 
■ ■ ■ there). Now notice that the integral in (6.46) becomes divergent when a > ^ for odd d 
and a > 1 for even d. In fact, even if the leading term in (6.46) is convergent, higher order 
terms in (6.44) will still eventually lead to divergent integrals (and will be more and more 
divergent). Physically such divergences reflect that higher order terms in the expansion of g 
are more irrelevant and thus are more sensitive to the short-distance structure near x = 0. 
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When (6.46) is convergent we find 



with 



1 ^ 1 



a < ^, d odd 



2r(|-») de^^n 
Now using odd for illustration, consider a > ^. Since the large R expansion of g 
assumes that jlz ^ 1, we should split the integral at x ~ i.e. 

dxg{xR)= j'^'' g{xR)+ j dxg{xR) . (6.49) 
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Since g{z) is well defined near z = 0, the first term in the above is convergent and is of order 
0{-^). Now plugging the large R expansion of g into the second term of (6.49) we find the 
previously divergent term becomes 

-— + ••• . (6.50) 



Similarly, the leading contribution of a term of the form z~"' with > 1 in the large z 
expansion of g{z) will always be of order 0{R^^) regardless of the value of n. We thus 
conclude that for a > | 

IdiR,a)--^ (6.51) 
llR 

where the numerical coefficient now receives the contribution from all orders of the large z 
expansion of g{z), and the validity of the above expression assumes that this infinite sum is 
convergent. When a < ^, then the leading dependence is given by (6.47) as potentially di- 
vergent higher order terms give only a contribution which is of order O(^) and is subleading 
compared to (6.47). 

Exactly the same can be said regarding (6.38) for even d except that the range becomes 
a > 1 and the leading order contribution becomes of order 0{R~^). 
We thus conclude that 

crm o(iR) , id~l)\\ eK bi&) , ji dodd 

^ ^ ^'^ ^ Id even 



30 



with b{a) given by (6.48), and 



{d-l)\\ eK «>^(oddd) 



= + . ' + Oie^) . (6.53) 



/^(w) «>l(evend) 



Note that equations (6.53) agree with the general expectation (2.14) of Sec. II. Here we find 
that when the IR dimension of the involved irrelevant operator is not too large one should 
instead have (6.52), which can be considered as predictions from the holographic duality. 
Finally when a = | for odd d and a = 1 for even d, similar argument leads to 

S (R) - 5^-) _,id-meKlH^ a=\ (odd d) 

[W « = l(evend) 

VII. SOME NUMERICAL STUDIES 

We now consider some specific holographic RG fiows, including some "realistic" ones 
which describe known field theory fiows at large and strong coupling. In these examples, 
the minimal surface cannot be found analytically and numerical calculations are needed. 
In some examples, the gravity solutions are only known numerically. In computing the 
entanglement entropy, one could to choose to evaluate (6.10) directly after obtaining the 
numerical solution for the minimal surface. The numerical integrations can sometimes be 
time-consuming. It is often more convenient to use (6.27) or (6.29) to obtain Sd directly 
from the coefficient Cd{R) in the asymptotic expansion of the minimal surface solution near 
the boundary. 

We will first consider d = 3, where we find that S^lR) is always non-negative and mono- 
tonic. This is no longer the case for d = 4. The non-monotonic examples include the GPPZ 
fiow [44] which describes the fiow of the Af = 4 SYM theory to a confining phase under a 
mass deformation and the so-called Coulomb branch fiow [45, 46] which describes the A/" = 4 
SYM theory in a Higgs phase. 

The two realistic examples below in Sec. VII A 1 and VII Bl have also been considered 
recently in [38]. 
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A. d = 3 



1. A realistic flow 

We first consider the holographic RG flow of [51] which describes the flow from M2 brane 
theory at UV to an A/" = 2 superconformal theory with an SU{'i) x f/(l) global symmetry 
at IR. The flow involves two scalar operators whose UV and IR dimensions are respectively 

A = (2,2), AiR= Q(l + v^),i(5 + v^)) . (7.1) 

Thus in (6.5) and (6.7), a = 1 and a = \{y/Yl - 1) ^ 1.562. Also = ^ ~ 1-299. The 
numerical results are presented in Fig. 3. Clearly, Sz{R) is positive, and monotonic. Note 
that for this flow, e in (6.30) is e ~ 0.3, and thus equation (6.37) should be a reasonable 
approximation, which is conflrmed by numerical results. Since a > ^, from (6.53) we expect 
S^{R) should fall off as 1/R for large R. This appears to flt the numerical data very well. 

2. A sharp domain wall: first order phase transition 

The earlier example is a rather shallow domain wall. Let us now consider a sharp domain 
wall. Unfortunately, there appears no realistic example of this type. We will thus play with 
toy examples by coming up with various monotonic functions f{z), a strategy which we will 
repeatedly use below. By scanning through various examples we flnd that even for -j^ ^ 1, 
S^lR) remains monotonic and approaches to the right asymptotic value. However, a new 
phenomenon arises when the domain wall becomes steep, as in the left plot of Fig. 4. In this 
case, the minimal surface action (6.10) can have more than one extrema for certain range of 
R, as indicated in the middle plot of Fig. 4. The minimal surface prescription instructs us 
to pick the one with the smallest area. While the entanglement entropy itself is continuous 
as a function of R, after taking derivative to obtain S^lR), we flnd there is a discontinuous 
jump at some value Rc as shown in the right plot of Fig. 4. In other words, there is a 
flrst-order "phase transition" in the entanglement entropy at R^. Despite the jump, S3 is 
still monotonic and approaches the right asymptotic value. 
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FIG. 3: Left: f{z) for the domain wall solution describing the flow of M2-brane theory to an 
IR fixed point preserving M = 2 supersymmetry. Middle: plot of S{R) — S\jy{R) where S'uv 
denotes that at the UV fixed point. The UV divergences cancel when taking the difference, but the 
resulting expression does not have a well-defined large R limit, with a linear /^-dependence. As in 
the case of a free massive scalar and Dirac field of Sec. V, the presence of such a linear term can be 
understood as a finite renormalization between the short distance cutoffs of the UV and IR fixed 
points, as discussed in Sec. II. Right: S3{R) for this flow. We normalize the value at UV to be 
1. The horizontal dashed line denotes the expected value (6.17) for the IR flxed point. The black 
line (lower line) is obtained from numerical calculation by using (6.27). For this flow, e in (6.30) is 
e w 0.3, and thus equation (6.37) should be a reasonable approximation, whose results are plotted 
using the red line (upper line). Note that the part linear in R in — 5finite as seen in the second plot 
is automatically eliminated when considering S^^R). 

3. "Gapped^' phases 

We now consider some examples where the IR is described by a "gapped" phase. We 
are interested in systems with Lorentz symmetry, and thus are forced to consider singular 
geometries discussed around (6.8). An example of gapped phase with a regular geometry 
is the AdS soliton [52]. But the corresponding boundary theory contains a circle direction. 
Our discussion of the divergence structure of the entanglement entropy in Sec. II thus does 
not apply, and one needs to modify the definition of Sd{R), which we will not pursue here. 
The holographic entanglement entropy for an AdS soliton was considered before in [53]. A 
first-order phase transition between minimal surfaces with disk and cylinder topology similar 
to our third example below (right plot of Fig. 5) was found there. 

Again due to lack of simple explicit examples, we will consider some toy examples by 
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FIG. 4: Left: a steep domain wall (toy example) with f(z) = 1 + 5100^^100 • Middle: plot of 
S{R) — SconiiR) where 5'conf(-R) denotes the entanglement entropy for the UV fixed point. The 
short-distance divergences cancel when taking the difi'erences. For the indicated range of R, the 
action (6.10) has three extrema, all of disk type. The entanglement entropy of the system is given 
by the smallest of them. There is a first-order "phase transition" at Rc = 4.4. Right: S3{R) has 
a discontinuous jump, which is indicated by the vertical green line. The dashed horizontal line is 
the expected asymptotic value for the IR fixed point. 



S,(R)/K Sm/K S,IR)/K 




R 



0.5 1.0 1,5 2.0 2.5 3.0 0.72 0.74 0.76 0.78 0.80 0.82 0.84 



FIG. 5: Left: SaiR) for f{z) = 1 + 2^. Middle: S^iR) for f{z) = 1 + 2^ which exhibits a "second- 
order phase transition" from minimal surface of disk topology (black curve) to cylindrical topology 
(red curve). Right: S3{R) for f{z) = (1 + z^)^, which exhibits a "first-order phase transition" 
between the surfaces of two topologies. The dashed curve corresponds to other extrema of the 
minimal surface action. There is a discontinuous jump in Ss{R) which is indicated by the green 
vertical line. 

postulating certain f{z). In Fig. 5, we consider three such examples which illustrate three 
possible scenarios. 

In the left plot we consider f{z) = 1 + z"^ . In this case since f{z — t- oo) ~ z"^ with an 
exponent 2, only minimal surface of disk topology is possible. There appears to be a unique 
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minimal surface solution for all i?, and S^lR) is smooth. 

The middle plot is for f{z) = 1 + z^, in which case, when R is sufficiently large, the 
minimal surface is of cylinder type (see discussion around (6.13) and Appendix C). This is 
indicated in the plot by the red curve. The black curve for smaller R values has minimal 
surface of disk type. S^lR) appears to be continuous at the transition point where the 
minimal surface changes topology. Thus the entanglement entropy has a "second-order 
phase transition." Note that approaching the transition point from the smaller R side, the 
end point of the disk-like surface approaches infinity, while from the larger R side, po 
for the cylinder-like surface approaches zero. The two branches meet at the critical surface 
discussed in Appendix C (see Fig. 6 for a cartoon). It may be interesting to understand the 
critical behavior of such a transition which we will leave for future study. 

The right plot of Fig. 5 is for f{z) = (1 + z^)^, which exhibits a "first-order phase 
transition." Again the red and black curves have minimal surface of cylinder-type and disk- 
type respectively. Now near the transition between the two topologies, the action (6.10) 
now has three extrema, with the non-minimal area extrema indicated in the plot by dotted 
lines. There is a discontinuous jump in S^lR) at the transition point. 

In all these examples, S^lR) appears to approach zero at large R, although our numerics 
cannot go to too large R. 

4- Non-monotonic f{z) 

Finally just for comparison, in Fig. 7 we show a toy model of / which is not monotonic, 
neither is the corresponding 1S3. This toy example illustrates nicely several important aspects 
of the properties of S^i^R): 

1. As mentioned earlier, the monotonicity of / is tied to the imposing of null energy 
condition, which can be interpreted as the bulk refiection of boundary unitarity. This 
indicates that the monotonicity of S^^R) is closely tied to unitarity of the boundary 
system. 

2. The location in R where the non-monotonicity of S^lR) occurs appears to roughly scale 
with z. From the IR/UV connection, z translats into a boundary length scale. That 
the non-monotonicity of f{z) at certain z directly translates into non-monotonicity of 
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FIG. 6: Cartoon of a "second-order phase transition" from a minimal surface of disk topology 
(black) to a minimal surface of cylinder topology (red) as we increase R. The surface (brown) at 
Rc is the critical surface discussed in Appendix C. 

f(Z) i'i(R)/K 



(S(R)-Suv(R))/K 
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FIG. 7: Left: A toy model of f{z) which is non-monotonic. Middle: The entanglement entropy 
itself is not sensitive to the non-monotonic feature of /. Right: the corresponding 53(i?) is sensitive 
and also becomes non-monotonic. 

Ss{R) at certain R which roughly scales with z indicates that SslR) is most sensitive 
to contributions from degrees of freedom around scale R. In contrast, as shown in the 
figure the entanglement entropy itself is rather featureless as a function of R. 
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B. d = 4 



1. Leigh- Strassler flow 

We first consider the flow from tlie A/" = 4 SYM tlieory (c? = 4) to an = 1 superconfor- 
mal flxed point witli SU{2) x f/(l) global symmetry [43, 54]. The flow involves two scalar 
operators whose UV and IR dimensions are respectively 

A = (2,3), AiR= (l + v^,3 + y7) . (7.2) 

Thus in (6.5) and (6.7), a = 1 and a = ^7 - 1 ^ 1.646. Also = ~ 1.12. The 
numerical results are given in Fig. 8. S^lR) is positive and monotonic. For this flow, e 
in (6.30) is e ~ 0.12, and thus equation (6.38) should be a good approximation, which is 
conflrmed by the numerical results. Since a > 1, from (6.53) we expect S4^{R) should fall 
off as l/R"^ for large R. Our numerical data for large R are not good enough to test this 
conclusively, but does not appear to directly contradict with it. 

f(z) 2S4(R)/K 

'■-[ 

1.0. ■ 

0.8 - 
0.6 - 

0.4- 
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2 4 6 8 

FIG. 8: Left: f{z) for the domain wall solution describing the Leigh-Strassler flow from AA = 4 
SYM theory to an IR fixed point preserving M = 1 supersymmetry. Right: S4{R) for this flow. 
The horizontal line denotes the expected value for the IR theory. The dotted line is obtained from 
direct numerical calculation. For this flow, e in (6.30) is e ~ 0.12, and thus equation (6.38) should 
be a good approximation, whose results are plotted using the solid line. 
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2. A sharp domain wall 

In the Leigh-Strassler flow discussed above, the central charges of the UV and IR fixed 
points are close. Our discussion in Sec. VI D indicates that in such a situation Sd{R) should 
be monotonic and positive. Now let us consider a toy example in which the domain wall is 
steep and the separation of central charges is large. As indicated in Fig. 9, Si{R) is neither 
monotonic nor positive definite. In contrast to the d = 3 example of Sec. VII A 2, there is 
no phase transition here. 

f(z) 



254(R)/K 




FIG. 9: Left: A steep domain wall (toy example) in d = 4 with f{z) = 1 + 515^^15 . Right: In 
contrast to the example of Sec. VII A 2 in d = 3, S4^{R) is smooth, but is neither monotonic nor 
positive definite. The dash line is the expected asymptotic value for the IR fixed point. 



3. Gapped phases I: GPPZ flow 

We now consider the GPPZ flow [44] which describes the flow of the A/" = 4 SYM theory 
to a confining theory under a mass deformation, which has UV dimension A = 3. For this 
flow the metric is known analytically with 

/W= (1 + 1^)' ■ (7.3) 

Note that as2;— J-oo, / j-oo. Low energy excitations of this system have a discrete 

spectrum with a finite mass gap (see e.g. [47]). For such a gapped phase we expect S/i{R) — )■ 
for large R. Fig. 10 gives the entanglement entropy and ^4 for this system. While the 
finite part of the entanglement entropy appears to grow linearly with R at large i?, 1S4 
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approaches zero from negative side. Similar to the f{z) = 1 + example in d = 3 discussed 
in Sec. V11A3, there is a "second order phase transition" from minimal surface of disk 
topology (black curve) to cylindrical topology (red curve), where S^^R) remains continuous 
(as far as our numerics could tell). Again as in the discussion of f{z) = 1 + z'^ the transition 
goes through the critical surface of Appendix C, as indicated in Fig. 6. 

2S4(R)/K 



(S(R)-Suv(R))/K 




FIG. 10: Left: plot of entanglement entropy S{R) — S\jy{R) for GPPZ flow where 5uv(^) 
denotes the entanglement entropy for the UV fixed point. The short-distance divergences cancel 
when taking the differences. Right: plot of 54(i?) for GPPZ flow, which is neither positive-definite 
nor monotonic. The system also exhibits a "second order phase transition" from minimal surface of 
disk topology (black curve) to cylindrical topology (red curve), where S4{R) remains continuous (as 
far as our numerics could tell). 



4- Gapped phases II: Coulomb branch flow 

The Coulomb branch flow describes the spontaneous breaking of the gauge symmetry (and 
also conformal symmetries) of the A/" = 4 SYM by a vacuum expectation value for an 
operator which is bilinear in the fundamental scalars [45, 46]. Thus here there is no operator 
deformation, and the end point describes a point on the Coulomb Branch. It involves a scalar 
field of UV dimension A = 3. The metric is also known analytically, although through an 
implicit function 

. = -(l-«)^«-. /=^^ (7.4) 
1^ 9v3 
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with II a mass scale characterizing the expectation values of the adjoint scalars of the A/" = 4 
theory. Near z = (boundary), 

/ = l + ^/iV + ---, (7.5) 

and near z = oo, 

f = If^'z' + ■■■ ■ (7.6) 

While the singularity z = oo lies at a finite proper distance away, but it takes a null geodesic 
an infinite time in t. The low energy excitations of the system include a massless Goldstone 
mode from spontaneous breaking of scaling symmetry, and a continuous spectrum above a 
finite mass gap (see e.g. [47]). Note that the single Goldstone mode will have a nontrivial 
0(1) contribution to ^4, but is not visible in the order 0{N'^) we are considering. So we 
expect S^{R) — )■ for large _R as in a gapped phase. With large z behavior given by (7.6), 
there can only be minimal surface with a disk topology. The numerical result is shown 
in Fig. 11, with a unique minimal surface solution for all R. However, S^lR) is neither 
monotonic nor positive-definite. 

2S4(R)/K 

1.0 
0.5 



FIG. 11: Si{R) for the Coulomb branch flow. 

C. Summary 

To summarize, we find that in d = 3 all examples (which are Lorentz invariant and 
satisfy the null energy condition) have a monotonic and non-negative S^lR), while this is 
no longer the case for S^^R) in ci = 4. We also find that the monotonicity of Ss{R) is tied 
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very closely to the null energy condition, which can be interpreted as the bulk reflection of 
boundary unitarity. Thus it is tempting to conjecture it is always monotonic in the vacuum 
of a Lorentz-invariant, unitary QFT. 

Note that although S^lR) is not monotonic, given the relation (4.11) and the recently 
proved a-theorem [32] for ci = 4, we still always have S4^{R — )■ 0) > 54(i? — )■ oo). 

Even if S^lR) turns out to be monotonically decreasing in the vacuum we do not expect it 
to remain so at a finite temperature/chemical potential. As discussed in Sec. II C at a finite 
temperature we expect Sd{R) to be proportional to the thermal entropy, i.e. proportional to 
the volume of the entangled region in the large R limit with a positive coefficient. Similarly 
for a system with a Fermi surface (i.e. at a chemical potential), we argued in Sec. Ill 
that Sd{R) should be proportional to the area of the entangled region at large R. While 
our dimensional analysis could not determine the prefactor, we expect it to be positive on 
physical ground. Clearly the corresponding Sd{R) which starts as a constant for small R, 
cannot not be monotonically decreasing in either situation. In the holographic context, 
the function f{z) in the spacetime metric for states of finite chemical potential or finite 
temperature does satisfy (6.4) (instead it decreases to zero at the horizon). 

We also observe first-order and second-order "phase" transitions in 5d(i?).^^ By first- 
order, we mean Sd{R) has a discontinuous jump, while in a second-order transition, Sd{R) 
is continuous, but not smooth. When the IR is described by a conformal fixed point (corre- 
sponding to a different AdS geometry), the phase transitions appear to be first-order and do 
not involve change of the topology of the minimal surface, as in a sharp domain wall in d = 3 
discussed in Sec. VII A 2. When the IR is described by a "gapped" phase corresponding to 
certain singular bulk IR geometries, the phase transitions (which can be first- or second- 
order) appear to involve change of topology of the minimal surface. Second-order transitions 
appear to be controlled by the critical surface discussed in Appendix C. Examples which 
exhibit a second-order transition include the GPPZ flow [44] which describes the flow of the 
A/" = 4 SYM theory to a conflning phase under a mass deformation. 

Such phase transitions are interesting, as they signal drastic changes of some underlying 
physics. In the example of a sharp domain wall, the system transitions very quickly (i.e. in 

Note that "phase" transitions in the holographic entanglement entropy have also been observed in various 
other contexts [53, 59-63], but all appear to be the first-order. 
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a very short range of scales) from the UV to IR regimes, thus it appears that the entan- 
glement entropy does not have "time" to respond to that quick change, which results in a 
discontinuity. Similarly the phase transitions in the case where the IR is a gapped phase 
likely signal that the system is opening a gap. Similar first-order phase transitions have been 
observed before in [60-62] and were interpreted as a confinement /deconfinement transition. 

VIII. CONCLUSIONS AND DISCUSSION 

In this paper we introduced a "renormalized entanglement entropy" (R) which ap- 
pears to capture only the "universal" part of the entanglement entropy. We illustrated the 
power of this construction by showing that the qualitative behavior of the entanglement 
entropy of a system with a Fermi surface could be obtained by simple dimensional analysis. 

We also showed that (in the vacuum) SP{R) has various nice features (as listed in 
the Introduction) which make it natural to interpret it as describing the RG flow of the 
entanglement entropy with distance scale. We were particularly interested in finding out 
whether it could be monotonic and provide a scale-dependent measure of the number of 
degrees of freedom of a quantum system. In = 3 our studies of free theories and holographic 
systems support the conjecture that S^lR) (given by (1.3)) is non-negative and monotonic 
for Lorentz-invariant, unitary QFTs. In ci = 4, this is no longer true. We find examples in 
which S^lR) (given by (1.4)) is neither monotonic nor positive-definite. 

Let us now mention some future questions: 

1. In our discussion (R) was constructed by using the differential operator in (1.1) 
to strip shorter-distance correlations (including cutoff dependence) away from the 
entanglement entropy S^^\R). In particular, to show its UV finiteness we had to 
examine the the specific structure of divergent terms in S^'^\R). It would be interesting 
to have an intrinsically finite way to define it. One possibility is that one may obtain 
it from certain limit of the mutual information (and their generalizations involving 
more than two regions). 

See e.g. [55, 56] for discussions mutual information in general quantum field theories and of possibility of 
a c-theorem from mutual information. 
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2. The monotonicity of (1.3) would imply that 

i.e. as a function R, S{R) is a concave function. This feels like a relation that can per- 
haps arise from a clever use of the strong subadditivity condition of the entanglement 
entropy, which is also responsible for the monotonicity of S2{R) in d = 2 [6]. It is also 
important to examine more examples for further confirmation or counterexamples. In 
the holographic context, for example, it would be interesting to generalize our discus- 
sion to gravity theories with higher derivatives as in [10, 34, 39-42, 57, 58] to test the 
robustness of the conjecture. 

3. Suppose Ssi^R) turns out to be monotonic, it is then rather curious this is no longer 
true for S^^R) in ci = 4. 

(a) One logical possibility is that S4^{R) is monotonic, but the holographic examples 
in which it is not are pathological. The example of Sec. VII B 2 is a toy model with 
a very sharp domain wall which is clearly artificial, while those of Sec. VII B 3 
and VII B 4 involve singular geometries. It would be important to completely 
settle this using more examples. 

(b) Now let us assume that S^lR) is indeed not always monotonic. Then one pos- 
sibility is that some modifications on the definition of 1S4 are needed in order to 
construct a central function out of the entanglement entropy on a sphere. An- 
other possibility is that the c-theorem for d = 4 may require physical information 
which is not encoded in the entanglement entropy for a spherical region. One 
hint which points to this thinking is that for = 4 or higher, the monotonicity 
of Sd{R) would involve more than two derivatives on S{R), e.g. for d = 4, it 
amounts to 

R:^dlS + R^dlS < ROrS , (8.2) 

while the strong subadditivity condition can only lead to the second derivative of 
S{R) when applied to infinitesimally separated regions. From this perspective, 
it appears likely Sd{R) will also not be monotonic for d > 5. It would be good 
to check this explicitly. It would also be interesting to understand whether non- 
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monotonic regions signal some interesting underlying physics at those distance 
scales. 



4. All examples of this paper have been on relativistic theories. It would be interesting 
to explore non-relativistic systems and also time-dependent systems. 

5. In this paper we find in holographic systems: 

(a) When the difference of the central charges of the UV and IR fixed points of a 
QFT is sufficiently small, then Sd is positive and monotonic for all d. 

(b) In the small R limit, the leading i?-dependence is given by 

Sd{R) = s^^""^ - + . . . (8.3) 

where a = c? — A for a source deformation and a = A for a vev deformation,^^ 
with A is the UV scaling dimension of the leading deformation operator (i.e. the 
one with the smallest a). A{a) is a positive number. 

(c) When the IR is described by a conformal fixed point, we expect in the large R 
limit, generically 

f,^. I H odd d 

SdiR) = + \ (8.4) 
I ^ + ■ ■ ■ even d 

with the sign of the coefficients undetermined. Bur for a < ^ for odd d and 
a < 1 for even d, where a = Air — d with Air the scaling dimension of the 
leading irrelevant operator in the IR, holographic systems predict that 

5,(i?) = sr + ^|p. + ---, R^oo (8.5) 

with B{a) a positive number. 

It would be desirable to derive these features more generally using conformal perturba- 
tion theory. It is also interesting to note that the asymptotic behavior (8.3) and (8.5) 
in fact coincide with those of the interpolating c-function of [32] (for d = 4).^^ Likely 
they have a common origin. 



This is for standard quantization. For alternative quantization one then have a = A for source deformation 
and a = d — A for vev deformation. 
^® We thank S. Pufu for this observation. 
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6. In both d = 3 and d = A we observed first-order and second-order "phase" transitions 
in Sd{R)- For a first-order transition, Sd{R) has a discontinuous jump, while in a 
second-order transition, Sd{R) is continuous, but not smooth. It would be interesting 
to understand whether such "phase transitions" are due to artifacts of the the large 
approximation one is working with. Even if they are, their presence should still reflect 
some underlying features of the system at finite A^. For example, as mentioned earlier, 
the transition in the GPPZ flow may signal the opening of a gap of the system. 
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Appendix A: Induced metric and extrinsic curvature for a scalable hypersurface 

Denote the spherical coordinates in M" as (r, 6a) where 6a,a = 1, ■ ■ - n — 1 denote all the 
angular variables. Then a scalable hypersurface S can be specified as 

r = RfiOa) , (Al) 

where f{Oa) is a smooth function of the angular variables only and R denotes the size. 
Clearly as we change R, the shape of the surface which is specified by / does not change. 
Plugging (Al) into the metric for M" we find that the full flat space metric for MJ^ can be 
written as 

ds'^ = ifdR + Rdff + R^fdn^ = fdR^ + 2RdRfdf + (A2) 
where dfl^ is the standard metric on a unit sphere and (is| is the induced metric on E 

dsl = R^ {fdQ^ + {dff) = habde^de^ . (A3) 

Note that since R^ appears in hab as an overall prefactor, the Christoffel symbol F built from 
hab is clearly independent of i?, and thus the intrinsic covariant derivative Da = de^+T on 
S is independent of R. From (A2) one can also read that the standard lapse function N is 
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independent of R, while the shift vector Na oc R. Thus the extrinsic curvature Kab which 
can be written as 

Kab = {dnKb - DaN, - D,Na) cc R . (A4) 

Thus we have estabhshed (2.5). Since in F of (2.1) all indices have to contracted and there 
are an even number of Kab, thus F can only depend on R through for n some non- 

negative integer. This then establishes (2.4).^^ Note that = 2 is a bit special as in this 
case, both the curvature K and h have only a single component, but the same conclusion 
applies. 

For illustration, let us also give an example of a surface which is not scalable. Consider 
in a curve specified by 



„2 , y _ ^2 
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(A5) 



whose shape clearly changes with R. For this curve there is an additional dimensional 
parameter b, and our previous discussion does not apply. 

Appendix B: Details of the numerical calculation of Ss{R) for a free massive scalar 



We start with the Hamiltonian for a free massive scalar in c? = 3 



(Bl) 



where 11 is canonical momentum for (p. In terms of the Fourier transform of (p and 11 in the 
angular 6 direction 




27r 



d0 



+ ic 



dO e''y{e,r), l>0 



Ho = ^ j de n(^, r), + in.i = ^ j de e'''u{e, r), / > o 

the Hamiltonian then can be written as H = X]«^-oo with 

f oo 



(B2) 
(B3) 



Uf + rdr 



(B4) 



For theories breaking parity it is possible to have terms with different powers from those indicated 
in (2.4) [13], but they have negative powers of R. From our discussion in the main text, they will 
not give rise to divergent terms. 
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Also note that 0/, 11; satisfy the canonical commutation relation 



i)i{r), Ui'ir')] = i5ii>6{r - r') . 



(B5) 



We discretize (B4) with a uniform lattice in the radial direction: 



Hj 



2a 



N 



2; L VJTT 



J' 



(B6) 



where a is the lattice spacing and r = ja and we introduced an IR cutoff Am = Na. The 
radius of the disk is taken to be: 



R= { n+^] a . 



The Hamiltonians (B6) can be written in a general form 



(B7) 



N 



N 



(B8) 



and we are interested in the entanglement entropy S{n, N) for the subset of degrees of 
freedom 0a,a; = l,2,---?T,<A^ when the system is in the vacuum. The problem was solved 
in [4] and the result can be written as follows. Decompose f2 = \/K as 



A B 
C 



(B9) 



where A has rank n and let 
1 



I3 = -B' A-'B 
^ 2 



-3- 



1 



/3' 



V2 



Then, S'(n, N) can be written in terms of the (A^ — n) x {N — n) matrix S as 



For our case: 



S{n,N) 



tr 



logfl - S) + 



IokS 



(BIO) 



(Bll) 



+ f + 



2 + 



J" 

J + 1/2 
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(B12) 



and the entanglement entropy for a spherical region of radius R can then be written as 

oo 

S{R)= Si{n,N) (B13) 

l=—oo 

where Si{n,N) is given by (Bll) with obtained from (B12) via (BIO) for each I. 

To get Si{n,N) we need the eigenvalues of S;. It can be checked that all square roots 
and inversions are well defined, and all the eigenvalues of lie between and 1. One can 
also show that the sum over I converges at fixed n, N. 

For m = 0, S{R) entropy should have the form 

SiR) = c^--s; + o(^,^) (B14) 
a \R AjrJ 



where S3 is known analytically [64] 



1 - 3C(3)\ 



S3 = — (^2 log 2 - J ^ 0.0638 . (B15) 

From our numerical calculations we get for m = 

ci = 0.4643821 S3 = 0.0635 ± 0.0004 (B16) 

The calculations were performed with the choice of IR cutoff = 200. To avoid the 
boundary effects from both the UV (small n) and IR (large n) cutoffs, we needed to restrict 
to the range 10 < n < 45. We require 10^^ absolute accuracy for the result of the entropy 
calculation in order to be able to extract SslR) with satisfactory precision. To achieve this 
we followed the analysis of [65] which shows that the finite volume corrections (accounting 
for a finite A^) go as A^^(2'+2). With our choice N = 200, they are thus negligible for / > 3. 
For / < 3 we determine the coefficient of the correction by doing the calculation for several 
choices of > 200 and extrapolate to infinite N. We refer the reader to [65] for details. 

At finite m, given the limited range of 10 < n < 45 (thus limited range of R) for which 
we can do calculations for a given m, we extend the range of mR by working with different 
choices of m in lattice units. Since S^lR) can only depend on mR in the continuum limit, 
data points from different choices of m should collapse into a single curve, and they indeed 
do. The results are shown in Fig. 1 of Sec. V. Choosing 1/m = 20a enabled us to cover 
the range 0.55 < mR < 2.25, while 1/m = 40a, 120a cover the range 0.275 < mR < 1.125 
and 0.09 < mR < 0.375 respectively. In Fig. 1 one can see that in the overlapping regimes 
the data points agree, justifying the continuum extrapolation of the lattice results, i.e. the 
S3{mR) determined here is independent of the discretization. 
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Appendix C: Cylinder-like solutions 



In this Appendix we discuss when a cyhnder-hke solution minimal surface could appear. 
We will show that for a bounded /, the minimal surface solution is always disk-like, while 
when f{z) ~ for large z, a cylinder-like solution can appear only when n > 2. 

Extremizing the minimal surface action (6.10) leads to the equation of motion 



Let us first consider a bounded /, i.e. with asymptotic behavior 

/(^) = /o + fiz-"' + ■■■, m>0, ^^oo. (C2) 

Then with a cylinder-like solution 

p{z) = po + Piz" + •••, a<0, z oo , (C3) 

in the equation of motion (CI), the leading term on the left hand side (LHS) is of order 0(1) 
while that on the right hand (RHS) is of order 0{z°'^'^), and thus a solution is not possible. 
Now let us consider an unbounded f{z) with the large z behavior 

/-I = az~" (1 + ■ ■ ■) , n>0, z^oo . (C4) 

With (C3), in the equation of motion (CI), the leading term on the left hand side is 

LHS = {d- 2)a2-" + {d - l)poPiaz"-'^ + ■■■ . (C5) 

The leading behavior on the RHS is given by 



RHS 



-|poPia(2 -n -2a)2-"-" a>l-f 

IpopMn + 2a- 2);z"-2 a< 1 - f (C6) 

~ z-^, /3>f + l a = l- f 



Comparing (C5) and (C6) we find that the leading order equation can only be satisfied for 
the case of the second line of (C6) with 

a = 2 - n . (C7) 
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Note that since a < 0, this requires n > 2. We also find that the equation of motion requires 



o = 2(rf-2)a 

Po{n-2){n + 2d-A) ' ^ ' 

We thus conclude that for n > 2, a cylinder-like solution could exist with the large z behavior 

given by 

PQ[n — 2)[n + 2d — 4) 

We notice that in (C9) the coefficient before the second term becomes singular when 
Po — ^ 0. This implies if there is a solution with po = it should comes with a larger power 
than (C7). Writing for large z, p{z) = p[z"' H — ■ with a' < 0, we again find that a solution 
is possible only for n > 2, for which the leading behavior of the solution is given by 



Equation (CIO) describes a solution in which the minimal surface just closes off at 2; = 00. 
Note that (CIO) can be considered as the po ~^ limit of (C9) in the following sense. 
Equation (C9) was derived assuming that the second term is much smaller than the first 
term, i.e. for 

n-2 1 2id-2)a 2 
' ^ plin-2)in + 2d-A)-'' ' 

When Pq ^ a, equation (CIO) is valid in the region 1 z < zq, while one has (C9) for 

z ^ zq. In the limit po — j- 0, 2:0 — )■ 00, and thus (CIO) becomes valid in the full large z 

region. 

The solution (CIO) can be considered as the critical surface lying at the boundary between 
the spaces of disk-like and cylinder-like solutions. As discussed in Sec. VII, indeed it appears 
to control the transition point of certain "second order phase transitions" from minimal 
surface of disk topology to cyhndrical topology. 



[1] L. Amico, R. Fazio, A. Osterloh and V. Vedral, Rev. Mod. Phys. 80, 517 (2008). 

[2] J. Eisert, M. Cramer, and M. Plenio, Rev. Mod. Phys. 82, 277 (2010). 

[3] L. Bombelli, R. K. Koul, J. Lee and R. D. Sorkin, Phys. Rev. D 34, 373 (1986). 

[4] M. Srednicki, Phys. Rev. Lett. 71, 666 (1993) [hep-th/9303048]. 



50 



[5] M. P. Hertzberg and F. Wilczek, Phys. Rev. Lett. 106, 050404 (2011) [arXiv: 1007.0993 [hep- 
th]]. 

[6] H. Casini and M. Huerta, "A Finite entanglement entropy and the c-theorem," Phys. Lett. B 
600, 142 (2004) [hcp-th/0405111]; 

H. Casini and M. Huerta, "A c-theorem for the entanglement entropy," J. Phys. A 40, 7031 
(2007) [cond-mat/0610375]. 

[7] A. B. Zamolodchikov, "Irreversibility of the Flux of the Renormalization Group in a 2D Field 
Theory," JETP Lett. 43, 730 (1986) [Pisma Zh. Eksp. Teor. Fiz. 43, 565 (1986)]. 

[8] S. Ryu and T. Takayanagi, "Aspects of Holographic Entanglement Entropy," JHEP 0608, 
045 (2006) [hep-th/0605073]. 

[9] J. L. Cardy, "Is There a c Theorem in Four-Dimensions?," Phys. Lett. B 215, 749 (1988). 
[10] R. C. Myers and A. Sinha, "Seeing a c-theorem with holography," Phys. Rev. D 82, 046006 
(2010) [arXiv:1006.1263 [hep-th]]; 

R. C. Myers and A. Sinha, "Holographic c-theorems in arbitrary dimensions," JHEP 1101, 
125 (2011) [arXiv: 101 1.5819 [hep-th]]. 
[11] D. L. Jafferis, I. R. Klebanov, S. S. Pufu and B. R. Safdi, "Towards the F-Theorem: N=2 
Field Theories on the Three-Sphere," JHEP 1106, 102 (2011) [arXiv:1103.1181 [hep-th]]. 

I. R. Klebanov, S. S. Pufu and B. R. Safdi, "F-Theorem without Supersymmetry," 
arXiv:1105.4598 [hep-th]. 

I. R. Klebanov, S. S. Pufu, S. Sachdev and B. R. Safdi, arXiv: 11 12.5342 [hep-th]. 
[12] D. Friedan and A. Konechny, Phys. Rev. Lett. 93, 030402 (2004) [hep-th/0312197]. 
[13] T. Grover, A. M. Turner and A. Vishwanath, "Entanglement Entropy of Gapped Phases 

and Topological Order in Three dimensions," Phys. Rev. B 84, 195120 (2011); 

[arXiv:1108.4038vl]. 

[14] M. Levin and X. G. Wen, "Detecting Topological Order in a Ground State Wave Function," 

Phys. Rev. Lett. 96, 110405 (2006) [arXiv:cond-mat/0510613]. 
[15] A. Kitaev and J. Preskill, "Topological entanglement entropy," Phys. Rev. Lett. 96, 110404 

(2006) [arXiv:hep-th/0510092]. 
[16] B. Swingle and T. Senthil, arXiv:1112.1069. 

[17] N. Ogawa, T. Takayanagi and T. Ugajin, arXiv:1111.1023 [hep-th]. 

[18] L. Huijse, S. Sachdev and B. Swingle, Phys. Rev. B 85, 035121 (2012); arXiv:1112.0573. 



51 



[19] E. Shaghoulian, arXiv:1112.2702 [hep-th]. 

[20] N. lizuka, S. Kachru, N. Kundu, P. Narayan, N. Sircar and S. P. Trivedi, arXiv: 1201. 4861 
[hep-th] . 

[21] S. S. Lee, Phys. Rev. B 80 165102 (2009); arXiv:0905.4532. 

[22] M. M. Wolf, Phys. Rev. Lett. 96 010404 (2006); [quant-ph/0503219]. 

[23] D. Gioev and I. Klich, Phys. Rev. Lett. 96 100503 (2006); [quant-ph/0504151]. 

[24] B. Swingle, Phys. Rev. Lett. 105, 050502 (2010) [arXiv:0908.1724]. 

[25] Yi Zhang, T. Grover, and A. Vishwanath, Phys. Rev. Lett. 107 067202 (2011); 
arXiv: 1102.0350. 

[26] C. G. Callan, Jr. and F. Wilczek, "On geometric entropy," Phys. Lett. B 333, 55 (1994) 
[hep-th/9401072]. 

[27] C. Holzhey, F. Larsen and F. Wilczek, "Geometric and renormalized entropy in conformal 
field theory," Nucl. Phys. B 424, 443 (1994) [hep-th/9403108]. 

[28] P. Calabrese and J. L. Cardy, "Entanglement entropy and quantum field theory," J. Stat. 
Mech. 0406, P06002 (2004) [hep-th/0405152]. 

[29] H. Casini and M. Huerta, "Entanglement entropy in free quantum field theory," J. Phys. A 
42, 504007 (2009) [arXiv:0905.2562 [hep-th]]. 

[30] A. CappcUi, D. Fricdan and J. I. Latorrc, Nucl. Phys. B 352, 616 (1991). 

[31] S. N. Solodukhin, "Entanglement entropy, conformal invariance and extrinsic geometry," Phys. 
Lett. B 665, 305 (2008) [arXiv:0802.3117 [hep-th]]. 

[32] Z. Komargodski and A. Schwimmer, arXiv: 1107.3987 [hep-th]; 
Z. Komargodski, arXiv: 11 12.4538 [hep-th]. 

[33] H. Casini, M. Huerta and R. C. Myers, JHEP 1105, 036 (2011) [arXiv:1102.0440 [hep-th]]. 

[34] L. Y. Hung, R. C. Myers and M. Smolkin, "On Holographic Entanglement Entropy and Higher 
Curvature Gravity," arXiv:1101.5813 [hep-th]. 

[35] M. Huerta, "Numerical determination of the entanglement entropy for free fields in the cylin- 
der," arXiv:1112.1277 [hep-th]. 

[36] S. Ryu and T. Takayanagi, "Holographic derivation of entanglement entropy from AdS/CFT," 
Phys. Rev. Lett. 96, 181602 (2006) [hep-th/0603001]. 

[37] T. Nishioka, S. Ryu and T. Takayanagi, J. Phys. A 42, 504008 (2009) [arXiv:0905.0932 [hep- 
th]]. 

52 



[38] T. Albash and C. V. Johnson, "Holographic Entanglement Entropy and Renormalization 

Group Flow," arXiv:1110.1074 [hep-th]. 
[39] J. de Boer, M. Kulaxizi and A. Parnachev, "Holographic Entanglement Entropy in Lovelock 

Gravities," arXiv: 1101. 5781 [hep-th]. 
[40] J. T. Liu, W. Sabra and Z. Zhao, "Holographic c-theorems and higher derivative gravity," 

arXiv:1012.3382 [hep-th]. 
[41] A. Sinha, "On higher derivative gravity, c-theorems and cosmology," arXiv: 1008.4315 [hep-th]. 
[42] M. F. Paulos, "Holographic phase space: c-functions and black holes as renormalization group 

flows," arXiv:1101.5993 [hcp-th]. 
[43] D. Z. Preedman, S. S. Gubser, K. Pilch and N. P. Warner, "Renormalization group flows 

from holography supersymmetry and a c theorem," Adv. Theor. Math. Phys. 3, 363 (1999) 

[hep-th/9904017]. 

[44] L. Girardello, M. Petrini, M. Porrati and A. Zaffaroni, Nucl. Phys. B 569, 451 (2000); 
[hep-th/9909047]. 

[45] D. Z. Preedman, S. S. Gubser, K. Pilch and N. P. Warner, JHEP 0007, 038 (2000) [hep- 
th/9906194]. 

[46] A. Brandhuber and K. Sfetsos, JHEP 9910, 013 (1999) [hep-th/9908116]. 

[47] M. Bianchi, D. Z. Frccdman and K. Skcndcris, JHEP 0108, 041 (2001) [hcp-th/0105276]. 

[48] L. -Y. Hung, R. C. Myers and M. Smolkin, "Some Calculable Contributions to Holographic 
Entanglement Entropy," JHEP 1108, 039 (2011) [arXiv:1105.6055 [hep-th]]. 

[49] C. R. Graham and E. Witten, Nucl. Phys. B 546, 52 (1999) [hep-th/9901021]. 

[50] A. Schwimmer and S. Theisen, "Entanglement Entropy, Trace Anomalies and Holography," 
Nucl. Phys. B 801, 1 (2008) [arXiv:0802.1017 [hep-th]]. 

[51] R. Corrado, K. Pilch and N. P. Warner, "An N=2 supersymmetric membrane flow," Nucl. 
Phys. B 629, 74 (2002) [hep-th/0107220]. 

[52] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998) [hep-th/9803131]. 

[53] A. Pakman and A. Parnachev, JHEP 0807, 097 (2008) [arXiv:0805.1891 [hep-th]]. 

[54] A. Khavaev, K. Pilch and N. P. Warner, "New vacua of gauged N=8 supergravity in five- 
dimensions," Phys. Lett. B 487, 14 (2000) [hep-th/9812035]. 

[55] H. Casini and M. Huerta, JHEP 0903, 048 (2009) [arXiv:0812.1773 [hep-th]]. 

[56] B. Swingle, arXiv: 1010.4038 [quant-ph]. 

53 



r] M. Fujita, arXiv:1112.5535 [hep-th]. 

i] N. Ogawa and T. Takayanagi, JHEP 1110, 147 (2011) [arXiv: 1107.4363 [hep-th]]. 
}] T. Hirata and T. Takayanagi, JHEP 0702, 042 (2007) [hcp-th/0608213]. 
3] T. Nishioka and T. Takayanagi, JHEP 0701, 090 (2007) [hcp-th/0611035]. 
L] I. R. Klebanov, D. Kutasov and A. Murugan, NucL Phys. B 796, 274 (2008) [arXiv:0709.2140 
[hep-th]]. 

2] I. Bah, A. Faraggi, L. A. Pando Zayas and C. A. Terrero-Escalante, arXiv:0710.5483 [hep-th]. 
3] M. Headrick, Phys. Rev. D 82, 126010 (2010) [arXiv: 1006.0047 [hcp-th]]. 
1] J. S. Dowker, "Entanglement entropy for odd spheres," [arXiv:1012.1548 [hcp-th]]. 
5] R. Lohmayer, H. Neuberger, A. Schwimmer and S. Theisen, Phys. Lett. B 685, 222 (2010) 
[arXiv:091 1.4283 [hep-lat]]. 
[66] R. C. Myers and A. Singh, "Comments on Holographic Entanglement Entropy and RG Flows," 
to appear. 



54 



